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ABSTRACT 



Aims. I investigate the shape of the electron cosmic ray spectrum in the range up to ~ 1000 keV, assuming that the acceleration 
process at the shock results in a power law in momentum, and that downstream of the shock the spectrum is affected by Coulomb 
interactions with background electrons only. 

Methods. In the non-relativistic regime one can analytically determine the energy of an electron starting with a certain energy, and 
use this result to produce an electron cosmic ray spectrum, modified by Coulomb losses. 

Results. An analytic expression for the electron spectrum is obtained that depends on the parameter n c t, which can be estimated from 
a similar parameter used to characterize the line spectra of supernova remnants. 

Conclusions. For the brightest supernova remnants n c t > 10" cirr 3 s, and most of the electrons accelerated to < 100 keV have 
lost their energy. Because of its high radio flux, Cas A is the most likely candidate for non-thermal bremsstrahlung. Although it has 
n c t ~ 2 x 10" cirr 3 s, one may expect to pick up non-thermal bremsstrahlung above 100 keV with current hard X-ray detectors. 

Key words. ISM: cosmic rays - acceleration of particles - supernova remnants - X-rays: ISM 



1. Introduction 

Over the last five years there has been considerable progress in 
our understanding of high energy cosmic ray acceleration by su- 
pernova remnants (SNRs). The spatial distribution and spectral 
shape of X-ray synchrotron emission from all historical SNRs 
and some other SNRs has made it clear t hat electron can be ac - 
celerated to energies up to ~ 100 TeV (Kov ama et al.L [19951) . 
that diffusive shock acceleration works very efficiently, w i th dif- 
fusion close to the Bohm limit dVinkl 120051: IVink et all 120061: 
iParizot et all l2006t IStage et all 120061) . and that magnetic fields 
are amplified by cosm i c ray streaming (IVi nk & Laming, 2003; 
iBerezhko et ail 120031: iBamba et all 120051: IVolk et all l2005l) . 
Moreover, the morphology of Tycho's SNR shows that at least 
in this remnant, but also probably in Kepler's SNR, cosmic ray 
acceleration is so efficient that it gives rise to enhance d com- 
pression ratios dWarren et all 120051: lEllison et all 120041) . In ad- 
dition, considerable progress has been made in the field of TeV 
astronomy, with many SNRs being e s tablished as sources of TeV 
y-rays (e.g. lAharonian etail 12004 120051: lAlbert et all 120071) . 
However, it has not yet been firmly established whether the y- 
ray emission is dominated by inverse Compton scattering from 
TeV electrons, or from neutral pion decay, caused by TeV ion 
cosmic rays. 

In contrast,little has happened concerning observations of 
low energy cosmic rays. This is unfortunate, since from a the- 
oretical perspective, the initial stages of the acceleration pro- 
cess, from thermal energies up to energies where Fermi accel- 
eration e fficiently operates, is co mplex and not well understood 
(see e.g. Ma lkov & Drurvil200lL for a review). In fact, the pro- 
cesses by which particles are injected into the Fermi acceler- 
ation process, are likely to be different for electrons and ions, 
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unlike Fermi acceleration itself. Our current understanding of 
the initial stages of acceleration is largely based on either in situ 
observations of inter-planetary shocks, and o n computer model- 
ing u s ing hybrid or particle in cell codes (e.g . lBvkov & UvarovL 
1999; SchmitzltiU l200l iLee et all l200l . There was some 
hope that hard X-ray observations of SNRs might provide obser- 
vational information on at least the electron component of low 
energy cos mic r ays dAsvarov et all, [T 990; Fav ata et all 119971 ; 
IVink etall 119971; IVink & Lamind.l2003h . However, it is gener- 
ally agreed that the non-thermal hard X-ray emission from SNRs 
such a s Cas A dThe et allll996tlAUen et allll997l:lFavata et all 
119971; IVink & Lamind. 120031: IRenaud et all l200 6bl). SN1006 
dAllen et all Il999t iKalemci et all 120061) Tycii o dAllen et all 
119991) . Kepler and RCW 86 dAllen et all 1 19991) is due to X- 
ray synchrotron from the highest energy electron cosmic rays, 
rather than from bremsstrahlung from low energy electrons. 
Nevertheless, the issues remains of interest as current and future 
hard X-ray detectors push the detection of hard X-ray emission 
to higher photon energies. At high enough energies, the steep- 
ening X-ray synchrotron spectrum is likely to be overtaken by 
non-thermal bremsstrahlung. 

In this paper I discuss the shape of the low energy cosmic ray 
electron spectra, given the importance of Coulomb-losses. These 
Coulomb losses in general alter the shape of the low electron 
cosmic ray spectrum as generated near the shock fron t. Note that 
a similar situation was investigated by Sarazin ( 1999) for higher 
energy electrons in clusters of galaxies. 

2. Coulomb loss affected electron spectra 

2. 1 . The electron energy as a function of time 

Supra-thermal electrons loose energy through various processes: 
bremsstrahlung losses, Coulomb losses (collisions with elec- 
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Fig. 1. The electron energy as a function of « e f (A ee = 30.9) 
for different values of energy at t = s (50 keV, 100 keV, 
200 keV...6400 keV). The solid lines show the exact, but im- 
plicit solutions to Eq.Q] i.e. Eq.|5] The dotted lines are the non- 
relativistic solutions to Eq.Q] i.e. Eq.|4] 
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Fig. 2. Low energy electron cosmic rays as expected from 
first order Fermi acceleration, including the effects of Coulomb 
losses (Eq. |5). The solid line gives the expected spectrum near 
the shock front, whereas the other lines show the models for dif- 
ferent values of n e f. 



trons/ions), and ionization lossesQ In the hot plasmas inside 
SNR shells Coulomb losses are likely to be the dominant source 
of energy losses of ele ctrons, in par ticular through electron- 
electron collisions (e.g. lHaugl 120041) . Electron-electron colli- 
sions are most efficient at low electron energies. For SNRs we 
are interested in the bremsstrahlung emission of supra-thermal 
electron with energies in the range of ~ 10- lOOOkeV. The prob- 
lem with this energy range is that neither the non-relativistic, 
nor the relativistic approximations for Coulomb losses are com- 
pletely valid. However, I will show that for most SNRs the non- 
relativistic approximation is sufficient. 

The energy loss rate fo r an el ectron with E » kT e is given 
by (e.g.Eiug, 2004;^uball2QQ2|): 



dE 2 , e + 1 

— = -Anrtm^c A ee n e , 

at p 



(1) 



with r e the classical electron radius, and e and p the kinetic en- 
ergy and momentum in dimensionless units, w R the electron den- 
sity, and A ee the Coulomb logarithm, given by (Huba, 2002): 



; ee = 30.9-ln[»>/ 2 (^p 
For the non-relativistic limit (e << 1) Eq. [Ogives: 



(2) 



dE 4nr 2 m 2 c 4 7.73 x 10~ 6 

-Aen-rie. = A ee ti e , (3) 



dt yjTm^E 



<E 



with energies, E in units of eV. 

The solution to this differential equation is: 

E(t) = (E 3 12 - 1.16 X 10 _5 ^ee«ef) 2/3 , 

with Eq the energy of the electron at t — 0. 



(4) 



1 In contrast, in the GeV-TeV range the dominant losses are due 
to synchrotron radiation and Coulomb l osses, limiting the maximum 
electron energy that can be obtained (e.g. lDrurv et al.L [l999). 



In fact, the correct equation, Eq. \T\ can also be solved, but 
only gives an implicit function for E(t): 



n e t ■ 



J(JK + l) 2 - 1 - arctan J(-^ + 1)2-1 

\ v mc 1 ' \ v mc- ' 



(5) 



+ l) 2 - 1 - arctan J(A + 1)2-1 



FigureQ]shows both the full solution and the non-relativistic 
approximation. As expected, the non-relativistic approximation 
fails at high energies. For energies < 400 keV the approximation 
is reasonable, with errors in the loss timescale of less than 30%. 



2.2. An analytic approximation to the non-thermal spectrum 

The advantage of an explicit function (Eq. |4]i is that it can be 
easily used to calculate the evolution of the electron spectrum, 
using the transformation N(Eo) — > N(E(t)), if an analytic ap- 
proximation is known. This can be done by inverting Eq.|4]and 
insert the result in N(E), taking into account that 



dEo 

He 



= (E V2 + 1.16 X 10-^ee«efr 1/J £ 



l/3pl/2 



(6) 



A simple example is when N(E) can be approximated by a 
power law N(Eo)dEo = KE r dEo. After some time t the spec- 
trum will be approximated by: 

N(E)dE = K(E 3 ' 2 + 1.16 x 10- 5 /i ee «efr (2r+1)/3 £ 1/2 d£- (7) 

A more interesting function in this context is the power 
law in momentum N(p) = Kp~ T , since the simplest version 
of first order Fermi acceleration predicts that this is the spec- 
trum of cosmic rays. Changing this into a function of energy 
(E = y/ p 2 c 2 + (mc 2 ) 2 - mc 2 ) gives (see lAsvarov et all 1 19901) : 

K r l 

N(E)dE = —r(E 2 + lErnc 2 )^ x (E + mc 2 )dE. (8) 
c 1 
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Inserting Eq.[4]and[6]in this expression gives: 

N(E)dE = -!-x (9) 
c 1 

((£ 3/2 + 1.16 X l(T 5 ,lee«ef) 4/3 + 

2mc 2 (E 3 ' 2 + 1.16 x l(T 5 ,l e e«e0 2/3 )~ <r+I)/2 x 
((£ 3/2 + 1.16 x 10" 5 ^ ee « e f) 2/3 + mc 2 )x 
(£ 3/2 + 1.16 x W- 5 A ee n e ty ll3 E l/2 dE. 

One of the important aspects of this solution is that it de- 
pends almost entirely on one parameter n e f, since A ee only 
weakly depends on the plasma density and temperature. This pa- 
rameter is familiar to anyone investigating X-ray emission from 
SNRs, as it is the key pa rameter for non-ionization equilibrium 
(NEI, see Liedahl, 1999) and for the equilibration of electrons 
and ions. In both cases the parameter is a measure for the num- 
ber of collisions with background electrons. 

Fig.|2]shows the electron distributions for values of n e t cov- 
ering the relevant range for SNRs, 10 9 - 10 12 cm" 3 s. Fig. Q] 
shows that for young SNRs the non-relativistic approximation 
is expected to be valid, as large errors only appear for n e f > 
5 x 10 11 cm~ 3 s. Moreover, the errors mostly affect the spectral 
shapes above 500 keV. 

3. Discussion and application to Cas A 

As explained above the shape of the non-thermal electron spec- 
trum below 1000 keV depends on the parameter n e f, which is 
the same parameter that governs the ionization state of the SNR. 
However, it depends on the details of the shock acceleration his- 
tory, whether one can safely assume that the average n e f of the 
hot plasma corresponds to the n e t governing the shape of the low 
energy electron cosmic rays. A close correspondence is likely if 
shock heating results always in a fixed ratio of heated and ac- 
celerated electrons at the shock front. Realistically, one should 
take into account the variation of n e t throughout the SNR, and 
the evolution with time. However, in practise a one or two val- 
ues for n e f characte rizes the overal l plasma characteristics rea- 
sonably well (e.g. lVinket~aTlll996l) . 

For young SNRs like Cas A, Tycho, Kepler and SN1006, the 
X-ray spectra are dominated by line emission from ejecta, shock 
heated by the reverse shock. This means that also n e t is deter- 
mined mostly by the plasma heated by the reverse shock, rather 
than by the forward shock heated plasma. It is often assumed that 
the reverse shock has a low magnetic field and therefore does 
not effi ciently accelerate cosmic rays. However, Hel der & Vinkl 
(2008) showed recently that the reverse shock of Cas A is re- 
sponsible for most of the X-ray synchrotron emission, this im- 
plies that the reverse shock is capable of accelerating particles 
as well.This is a confir mation of recent theoretical calculations 
bv lEllison et al. I (120051) . Nevertheless, it is not quite clear what 
the ratio is of the electrons accelerated by the forward shock, 
as opposed to the reverse shock. Moreover, the initial accelera- 
tion process for electrons may be different: The reverse shocked 
plasma is dominated by metals. In that case, a sizable fraction 
of the accelerated electrons may originate from ionization of 
(mildly) accelerated ions. 

Apart from the plasma « e f there is another reason why the 
n e t for the accelerated electrons is expected to be high in Cas A: 
Most of the electrons seem to have been accelerated in the early 
life of the SNR. The best evidence for this is the ~ 1 %/yr decline 



in its radio flux, which is best explained by adia batic expansion 
of the relativistic electron gas (Shklovsky, 1968). Cas A may be 
a special case since it is evolving into a red super giant wind, 
which has a density falling of with radius as r" 2 . This means 
that the particle flux entering the forward shock at any radius is 
more or less constant, whereas the shock velocity is dropping 
with radius. As a result particle acceleration was more efficient 
early on, suggesting a high n e t. 

Unfortunately, this also means that it is quite unlikely that 
the non-thermal X-ray emission of Cas A below ~ 100 keV is 
dominated by non-thermal bremsstrahlung, either fro m the for- 
ward s h ock, or from internal sh ocks, as assumed by (Laming, 
12001 alibi IVink & LamingLl2003l) . The reason is that all the emis- 
sion must then come from a limited fraction of the total shock 
heated plasma for which n e t < 10 10 cm~ 3 s (see[Willingale et al., 
120021; lYanget all 120071 for the spatial distribution of n e t ). If I 
simply assume a linear relation between age and n e t, less than 
5% of the plasma has n e t < 10 10 cm" 3 s. The lower hybrid wave 
model propagated by Laming] d2001al) remains an attractive op- 
tion for the injection of electrons into the Fe rmi acceleration pro- 
cess itself (see also Ghavamia n et al.L [20071) . 

Nevertheless, Cas A remains one of the most important tar- 
gets for searching for non-thermal bremsstrahlung, but only 
above ~ 100 keV, where the spectrum is less affected by 
Coulomb losses. The reason is that Cas A is by far the brightest 
Galactic SNR in the radio. This implies the ample presence of 
relativistic electrons, in com bination with a relatively high mag- 
netic f ield of 100 - 500 uG dShklovskvi 119681: 1 Vink & Lam ing. 
120031: iBerezhko & vEu l200l) . This is illustrated in Fig. [3] 
which suggests that with INTEGRAL-ISGRI we may be close 
to a detection, or even that part of the flux density around 
100 keV is caused by non-thermal b remsstrahlung (c.f. the 
CGRO-OSSE detection. lThe et all 1 19961) No d eep Suzaku Hard 
X-ray Detector (HXD) (Takahashi et all 120071) of Cas A exist, 
but a deep observation with this sensitive instrument, may also 
be able to detect emission above 100 keV. 

The normalization of the spectrum in Fig. [3] is determined 
from the radio flux, and assuming a mean magnetic field of 
~ 300 fiG. The bremsstrahlung also depends on the emission 
measure: EM = J E,H,H e Z 2 (f V, with n e being, here, the den- 
sity of the accelerated electrons. For Cas A the accelerated elec- 
trons are present in both the forward shock heated plasma, with 
«h ~ 10 cm" 3 , and in the reverse shock heated plasma, which 
consists of pure metals, p robably dominat ed by almost com- 
pletely ionized oxygen (e.g. lVink et all[l996l) . The emitting vol- 
ume of the forward shock heated plasma is Vf ~ 1.2 x 10 57 cm 3 , 
based on a shock radius of 2.6 pc and a contact discontinuity at 
~ 2 p c. The reverse shock radius is ~ 1.9 pc dHelder & Vinkl 
2008), thus V r ~ 1.4 X 10 56 cm 3 . Conservatively assuming that 
there is — 1 M Q of ionized oxygen dVink et all 1 19961) and that the 
other metals contribute little to the bremsstrahlung, one finds for 
the shocked ejecta 2,n,Z 2 w 35 cm" 3 . For the volume averaged 
normalization one finds (V/nH + V,-2/n;Z 2 )/(V/ + V r ) « 13 cm -3 . 
Note that both the bremsstrahlung and synchrotron luminosity 
scale with the emitting volume. It is not unreasonable to assume 
that the synchrotron emitting volume is equal to the volume of 
the shock heated plasma. However, it is uncertain how the ac- 
celerated electron density and magnetic field is distributed over 
this volume. So the normalization of < n,Z 2 >= 10 cm" 3 
used in Fig.[3]should be considered as a reasonable guess. For a 
constant non-thermal bremsstrahlung emissivity a higher value 
of S < «,Z 2 > implies a higher magnetic field. 
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Fig. 3. Left: Spectral energy distribution of Cas A, as measured with BeppoSAX-PDS (IVink et all 1200 ll) and INTEGRAL-IBIS 
dRenaud et alll2006bl) . assuming a magnetic field of ~ 300 //G, a spectral energy index of -2.56 (corresponding to a radio-spectral 
index of a — -0.78), and a backgrou nd plasma with Z, < «,Z 2 >= 10 cm 3 . The bremsstrahlung spectra were calculated using the 
analytic cross sections of lHaugl(ll997l) . The model spectra are shown for n e t — 0, 2 x 10 9 , 2 x 10 10 , 2 x 10 11 and 2 x 10 12 cnT 3 s (from 
top to bottom). Right: Spectral energy distribution of Tycho's SNR, as measured by the BeppoSAX-PDS. The theoretical curves are 
for B = 10 juG, a = -0.6, and E; < n,Z 2 >= 4 cm -3 . The different curves are for n e t = 0, 1 x 10 10 , 1 x 10 11 cm~ 3 s. 



In contrast, for Tycho's SNR, due to its much lower ra- 
dio luminosity, one only expects a detection of the non-thermal 
bremsstrahlung, if the magnetic field turns out to be unexpect- 
edly low: B < 10 //G This is roughly the value of the com- 
pressed average ISM magnetic field, whereas estim ates for the 
magn etic field in Tycho indicates B ~ 300 fiG dVolk et all 
l2005h . 

Finally, one could also consider SNRs in which the den- 
sities and n e t values are modest, and for which also the mag- 
netic fields are lo w; for example SN1006 (n e ? = 2 x 10 9 cm" 3 s 
IVink et all |2003l) and the Northeast of RCW 86 (n e t ~ 5 x 
10 9 cm s IVink et a 1., 2006). The advantage is that although 
the bremsstrahlung may be weaker, the low n e t ensures non- 
thermal bremsstrahlung at lower photon energies. Unfortunately, 
both these SNRs have also X-ray synchrotron radiation, which 
may be hard to distinguish from non-thermal bremsstrahlung. 
However, with imaging spectroscopy above 10 keV one may 
isolate regions where non-thermal bremsstrahlung is important 
and X-ray synchrotron radiation absent. Moreover, high spectral 
resolution spectra may reveal departures from Maxwellian elec- 
tron distributions. This type of research will have to wait till the 
emergence of high resolu tion, high thr oughput, X-ray imaging 
spectroscopy with XEUS dParmar et alll2006l) . 



4. Summary and conclusions 

I have derived an analytic expression for the low energy electron 
cosmic ray spectrum, assuming that the cosmic ray spectrum, as 
produced by the diffusive acceleration process, is a power law in 
momentum, and is only affected by Coulomb losses downstream 
of the shock. As may be expected, the parameter governing the 
Coulomb losses depends on « e f, the product of (background) 
electron density and the age of the accelerated electron popu- 
lation. This parameter is similar to the ionization time scale of 
the hot, shock heated plasma, and can be easily estimated from 
the line spectra of SNRs. 

For the SNRs with high densities, and large n e t, one does 
not expect a large population of accelerated electrons below 
100 keV. This is in particular the case for the very bright ra- 



dio source Cas A. Nevertheless, Cas A remains a good candi- 
date for searching for non-thermal bremsstrahlung, because of 
its high ion density and because its radio brightness indicates 
a large amount of accelerated electrons. Based on reasonable, 
but somewhat uncertain, assumptions, the results presented here 
imply that only above ~ 100 keV one expects to pick up the 
non-thermal bremsstrahlung component from Cas A. This may 
be accomplished with future deep Suzaku HXD observations, or 
additional deep INTEGRAL observations. 
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